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ABSTRACT 

j> . We investigate the deposition of energy due to the annihilations of neutrinos 

ON 



(N 



and antineutrinos on the rotation axis of rotating neutron and quark stars, 
respectively. The source of the neutrinos is assumed to be a neutrino-cooled 
accretion disk around the compact object. Under the assumption of the sepa- 



O 
O 



rability of the neutrino null geodesic equation of motion we obtain the general 
relativistic expression of the energy deposition rate for arbitrary stationary 
and axisymmetric space-times. The neutrino trajectories are obtained by us- 
ing a ray tracing algorithm, based on numerically solving the Hamilton- Jacobi 
equation for neutrinos by reversing the proper time evolution. We obtain the 
energy deposition rates for several classes of rotating neutron stars, described 
by different equations of state of the neutron matter, and for quark stars, 
described by the MIT bag model equation of state and in the CFL (Color- 
Flavor-Locked) phase, respectively. The electron-positron energy deposition 
rate on the rotation axis of rotating neutron and quark stars is studied for 
two accretion disk models (isothermal disk and accretion disk in thcrmody- 
namical equilibrium). Rotation and general relativistic effects modify the to- 
tal annihilation rate of the neutrino-antineutrino pairs on the rotation axis of 
compact stellar, as measured by an observer at infinity. The differences in the 
equations of state for neutron and quark matter also have important effects on 
the spatial distribution of the energy deposition rate by neutrino-antineutrino 
annihilation. 
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1 INTRODUCTION 
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The study of the electron energy deposition rate from the v + v — > e + + e~ neutrino an- 
nihilation reaction was initiated by Cooperstein et al. (1986), Cooperstein et al. (1987), 
and Goodman et al. (1987), respectively, and this process has been intensively investigated 
in the physical and astrophysical literature. Neutrino-antineutrino annihilation into elec- 
trons and pos itrons can deposit mo re than 10 51 ergs above the neutrino-sphere of a type 
II supernova (IGoodman et al. 19871 ). For a full understanding of the effects of the neu- 
trino annihil ation in strong gravitationa l fields, general relativistic effects must be taken 
into account ( jSalmonson fc Wilson 19991 ) . In a Schwarzschild geometry, the efficiency of the 
v + v — > e + + e~ process is enhanced over the Newtonian values up to a factor of more than 
4, in the regi me applicable to Type II sup ernovae, and by up to a factor of 30 for collapsing 
neutron stars (jSalmonson fc Wilson 19991 ). The neutrino pair annihilation rate into electron 
pairs between two neutron stars in a binary system was calculated by Salmonson & Wilson 
(2001). 

The gravitational effects on neutrino pair annihilation near the neutrinosphere and 
around the thin accretion disk were considered in Asano & Fukuyama (2000), by assuming 
that the accretion disk is isothermal, and that the gravitational field is dominated by the 
Schwarzschild black hole. General relativistic effects were studied only near the rotation 
axis. Using idealized models of the accretion disk, Asano & Fukuyama (2001) investigated 
the relativistic effects on the energy deposition rate via neutrino pair annihilation near the 
rotation axis of a Kerr black hole, by assuming that the neutrinos are emitted from the 
accretion disk. The Kerr parameter, a, affects not only the behavior of the neutrinos, but 
also the inner radius of the accretion disk. When the deposition energy is mainly contributed 
by the neutrinos coming from the central part, the redshift effect becomes dominant as a 
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becomes large, and the energy deposition rate is reduced. On the other hand, for a small 
a, the bending effect becomes dominant and makes the ene rgy increase by factor of 2, as 
compared with the case that neglects the relativistic effects (lAsano fc Fukuyama 200ll ). 

General relativity and rotation cause important differences in the spa tial distribution o f 
the energy deposition rate by neutrino v and antineutrino //-annihilation (IBirkl et al. 20071 ). 
The energy-momentum deposition rate (MDR) from the v — v collisions above a rotating 
black hole/thin accretion disk system was calculated by Miller et al. (2003), by imaging the 
accretion disk at a specified observer using the full geodesic equations, and calculating the 
cumulative MDR from the scattering of all pairs of neutrinos and antineutrinos arriving at 
the observer. The dominant contribution to the MDR comes from near the surface of the 
disk with a tilt of approximately 7r/4 in the direction of the disk's rotation. The MDR at 
large radii is directed outward in a conic section centered around the symmetry axis and is 
larger by a factor of 10-20 than the on-axis values. There is also a linear dependence of the 
MDR on the black hole angular momentum. The deposition of energy and momentum, due 
to the annihilation of neutrinos v and antineutrinos v in the vicinity of steady, axisymmetric 
accretion tori around stellar- mass black holes was investigated in Birkl et al. (2007). The 
influence of general relativistic effects were analyzed in combination with different neutri- 
nosphere properties and spatial distribution of the energy deposition rate. Assuming axial 
symmetry, the annihilation rate 4-vector was numerically computed. As compared to New- 
tonian calculations, general relativistic effects increase the total annihilation rate measured 
by an observer at infinity by a factor of two when the neutrinosphere is a thin disk, but the 
increase is only 25% for toroidal and spherical neutrinospheres. Thin disk models yield the 
highest energy deposition rates for neutrino-antineutrino annihilation, and spherical neutri- 
nospheres the lowest ones, independently of whether general relativistic effects are included 
or not. 

The study of the structure of neutron star disks based on the two-region (i.e., inner and 
outer) disk scenario was performed by Zhang & Dai (2009), who calculated the neutrino an- 
nihilation luminosity from the disk in various cases. As compared with the black hole disk, 
the neutrino annihilation luminosity above the neutron star disk is higher. The neutron star 
disk with the advection-dominated inner disk c ould produce the hi ghest neutrino luminosity, 
while the disk with an outflow has the lowest ( Jzhang fc Dai 2OO9I ) . A detailed general rela- 
tivistic calculation of the neutrino path for a general metric describing a rotating star was 
studied in Mallick & Majumder (2009). The minimum photosphere radius was calculated 
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for stars with two different equations of state, each rotating with two different velocities. 
The results show that the minimum photosphere radius for the hadronic star is greater than 
for the quark star, and that the minimum photosphere radius increases as the rotational 
velocity of the star decreases. The minimum photosphere radius along the polar plane is 
larger than that along the equatorial plane. The estimate of the energy deposition rate for 
neutrino pair annihilation for the neutrinos coming from the equatorial plane of a rotating 
neutron star was calculated along the rotation axis in Bhattacharyya et al. (2009), by us- 
ing the Cook-Shapiro- Teukolsky metric. The neutrino trajectories, and hence the neutrino 
emitted from the disk, are affected by the redshift due to the disk rotation and gravitation. 
The energy deposition rate is very sensitive to the value of the temperature, and its varia- 
tion along the disk. The rotation of the star has a negative effect on the energy deposition 
rate, it decreases with increase in rotational velocity. A phase transition during a super- 
nova explosion ca n induce stellar collapse and result in large amplitude stellar oscillations 
(IChan et al. 20091 ) . Extremely intense pulsating neutrino fluxes, with submillisecond period 
and with neutrino energy (greater than 30 MeV), can be emitted because the oscillations 
of the temperature and density are out of phase almost 180°. The dynamical evolution of 
a phase-transition-induced collapse of a neutron star to a hybrid star, which consists of a 
mixture of hadronic matter and strange quark matter, was studied in Cheng et al. (2009). It 
was found that both the temperature and the density at the neutrinosphere are oscillating 
with acoustic frequency. Consequently, extremely intense, pulsating neutrino/antineutrino 
fluxes will be emitted periodically. Since the energy and density of neutrinos at the peaks 
of the pulsating fluxes are much higher than the non-oscillating case, the electron/positron 
pair creation rate can be enhanced significantly. Some mass layers on the stellar surface can 
be ejected, by absorbing energy of neutrinos and pairs. These mass ejecta can be further ac- 
celerated to relativistic speeds by absorbing electron/positron pairs, created by the neutrino 
and antineutrino annihilation outside the stellar surface. 

It is the purpose of the present paper to consider a comparative systematic study of 
the neutrino-antineutrino annihilation process on the rotation axis of rotating neutron 
and strange stars, respectively, and to obtain the basic physical parameters characteriz- 
ing this process (the electron-positron energy deposition rate per unit volume and unit 
time), by taking into account the full general relativistic corrections. The source of the 
neutrinos is assumed to be th e accretion disk that can be formed around neutron and 
quark stars (IZhang fc Dai 20091 ). Neutrino-cooled accretion disks around stellar mass black 
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holes or neutron sta rs are plausible candidates for the central engine of gamma-ray bursts 
( jZhang fc Dai 20091 ) . As a first step in obtaining the electron-positron energy deposition 
rate on the rotation axis of different types of neutron and quark stars we generalize the 
relativistic description of the neutrino-antineutrino annihilation process to the case of arbi- 
trary stationary and axisymmetric geometries. By assuming that the neutrinos move on null 
geodesies, the main kinematic parameter of the collision process (the scalar product of the 
four momenta of the colliding vv pairs) are derived under the assumption that, similarly to 
the case of the Kerr black holes, the geometric variable in the geodesic equation of motion 
for the neutrinos can be separated. The domain of validity of this assumption is carefully 
investigated, and it is shown that in the case of the low mass stellar models with masses 
around M = 1.8M , the error is under 0.5% for angular velocities of the order of 2 x 10 3 s -1 , 
while it is more the 1% for angular velocities higher than 4 x 10 3 s _1 . The electron energy 
deposition rate is calculated for a general axisymmetric metric under the separability as- 
sumption. The neutrino trajectories are studied by using a ray tracing algorithm, based on 
numerically solving the Hamilton- Jacobi equation for neutrinos by reversing the proper time 
evolution. The electron-positron energy deposition rate on the rotation axis of rotating neu- 
tron and quark stars is studied for two accretion disk models (isothermal disk and accretion 
disk in thermodynamical equilibrium). In all these cases the electron energy deposition rate 
is obtained for different equations of state of the neutron and quark matter, as well as for 
different physical parameters of the stars. 

In order to compute the electron-positron energy deposition rate, the metric outside the 
rotating general relativistic stars must be determined. In the present study we obtain the 
equilibrium co nfigurations of the rotating neu tron and quark stars by usin g the RNS code, a s 
introduced in ( jStergioulas fc Friedman 19951 ). and discussed in detail in (ISterrioulas 20031). 
This code was used for the study of different models of rot ating neutron stars in ([Nozawa et al. 19981 ) 
and for the study of the rapidly rotating strange stars ( jStergioulas et al. 19991 ) . The soft- 
ware provides the metric potentials for various types of compact rotating general relativistic 
objects, which can be used to obtain the electron-positron energy deposition rate on the 
rotation axis of rotating neutron and quark stars. 

The present paper is organized as follows. In Section |2] we present the basic formalism 
for the calculation of the kinematic parameters (four-momenta) of the neutrino-antineutrino 
annihilation. The general relativistic energy deposition rate is obtained in Section |3j The 
equations of state of dense neutron and quark matter used in the present study are presented 
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in Section HJ In Section [5] we obtain the electron-positron energy deposition rates on the 
rotation axis of the considered classes of neutron and quark stars by assuming that the disk 
is isothermal. The case of the disk in thermodynamic equilibrium is discussed in Section |6j 
We discuss and conclude our results in Section [7J 



2 FOUR-MOMENTA OF MASSLESS PARTICLES AT THE ROTATIONAL 
AXIS IN STATIONARY AND AXISYMMETRIC SPACETIMES 

Let us consider in the coordinate system (t, r,9,(f>) an arbitrary stationary and axially sym- 
metric spacetime, given by the metric 

ds 2 = c 2 g tt dt 2 + cg t(j) dtd(j) + g„-dr 2 + geed9 2 + g^dcf) 2 , (1) 

and adapted to the symmetries of the geometry, with the time- and space-like Killing vec- 
tor fields (d/dt) and (d/d<fi), respectively. In the adapted coordinate system the metric 
coefficients depend only on the coordinates r and 9, measuring the curvature, isotropic, 
quasi-isotropic etc. radial distances, and the poloidal angle, respectively. 

For massless particles propagating along null geodesies in the spacetime (PQ), the geodesic 
equations containing the time-derivatives r = g rr p r and 9 = g e9 pe can be written as 
.rrji , „eej2 _ 

9t<t> ~ 9tt94>4> 

with the constants of motion p t = — Uq/c and p^ = L. Here Pt-.Pr-.Pe and p^ are the covari- 
ant components of the 4-momentum p of the light like particles. By applying the formal 
representations 

9 rr 9ee = 9ee / g rr = f (r) , (3) 
g2{r) -a 2 sin 2 9, (4) 



rr 2 , ee 2 _ 9^o/c) 2 + 2g t(j) (uj /c)L + g tt L 2 
9 P r i 9 Ve — 2 ' \ > 



9 2 t(i> ~ 9u9h> f( r 

gt^gee g(r)h(r,9) 



a — 



9t4> ~ 9u9h f{r) 
gttgee h 2 (r,9) 1 



(5) 



(6) 



g 2 ^-gttgu f( r ) sin2 #' 

Eq. (J2J) can be written as 

^H(^- fti )H>*-^) 2 - < 7) 

In Eqs. we have introduced a constant a, and the functions / = f(r), g — g(r) and 

h = h(r,9), and we have assumed that h depends on both r and 9, but / and g are only 
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functions of the radial coordinate. This is a rather strong assumption, since these functions 
are constructed from the metric coefficients via Eqs. ©-(El), which for a general stationary 
and axially symmetric spacetime depend on both r and 9, 
1 



g± 
h 2 



9tt 



9ee 



9tt9ee 



+ 



1 



(8) 



(9) 



g rr \9t4,- 9tt94><t> sin J 
We will examine the validity of these assumptions for the neutron and quark star spacetimes 

under our study in Appendix A. In the following we assume that g(r) ~ g+{r, 0) is a good 

approximation for the function g. By virtue of the pure radial dependence of / and h, the 

geodesic equation ((Tj) for L = can be separated with respect to the coordinates r and 9: 



Pr 



g+ u 


f-i 


f c, 




( wo 


sin 9 I 


a — 


V c 





J 



(10) 

(11) 



where J is a separation constant. In the following we use the positive branch g + of Eq. (JHJ), 
which gives r 2 + a 2 for the Kerr metric, where a becomes the rotational parameter, and leads 
to consistent identities in the decompositions in Eqs. (jlj) and (j3]). If we insert the expressions 
for p r and pg given by Eqs. fflO|) and fflTl) into the scalar product of the 4-momenta of a 
scattering vv pair for L = we obtain 



Pu Pu 



g 



c c 



+ g PvrPvr + g PvdPve 



-g +g 



c c 



+g 




ee LO 0v U 0P C 2 J 2 



a 2 sin 2 9. 



By introducing the variable p u , and the horizontal collision angle 9 V 



cJu . a V7 

, sin 9 V = p v , 

uou g+ 



we can cast the result in the form 

Pv Pv 



Cf 



g tt + g ee9± cosdv cose _ +g 6 



\ f 



sin 2 9 V -a 2 sin 2 9, 



g\ 



\ f 



sin 2 9 D - a? sin 2 9 



In the following we restrict our study to the neutrino-anti neutrino pairs that are scattering 
at the rotational axis. For 9 = the scalar product reduces to the expression 



P v Pu 



^ 2_ 



g lt + g 



ee g 



9a 



—t- cos 9 V cos 9 P + sin 9 V sin 9 P 
J J j 
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— g -jj- [1 — cos 9 V cos Op — sin 9 V sin fl c 



c- " / 

where we have used the relation g** = —g^/igf^ — 9tt9<t><t>) an d Eq. (jl]), evaluated at 9 = 0. 
By introducing the proper energy of the neutrinos defined as 



we obtain 



9t(j> ~ 9tt94><t> 



0„ 



6»=0 



\ 



9 9e (r,0) 



9l(r,0) 
f(r) 



[1 — cos 6 V cos 9 P — sin 9 U sin £ 



or, in the general case 



Pv Pv 



[1 — COS6 1 !, cos 6* p — sin^ sin^p cos(<^ !/ — <^p)], 



(12) 



(13) 



where ip u and ip„ are the azimuthal collision angles for the vv pair, as measured at the 
rotational axis. 



This equation is the same as the one derived for the Kerr spacetime (lAsano &: Fukuyama 20011 1 
Provided that neglecting the angular dependence of the function g + , and the substitution 
of the separated geodesic equations ffTUl) and ffTTj) into the scalar product does not increase 
the error propagation significantly, the expression (fl3|) is a good approximation for any 
stationary and axially symmetric spacetime. We can therefore use the same procedure to 
calculate the neutrino-antineutrino annihilation energy deposition rate at the rotational axis 
of neutron and quarks stars as the one used for the black holes. 



3 ENERGY DEPOSITION RATE FROM vv PAIR ANNIHILATION 
ALONG THE ROTATIONAL AXIS 

We consider a geometrically thin accretion disk around a rotating compact object, with 
an inner edge at the curvature radius R in , and a cutoff radius at Rout- We assume that 
the mass-energy of the disk produces only a negligible effect on the spacetime geometry 
of the central object. Thus any massless particle can propagate along null geodesies in the 
vicinity of the disk-compact object system, with the global geometry determined by the 
central massive object. We also suppose that the disk is the primary source of neutrinos in 
the system, and we neglect any interaction between the disk and the neutrino/anti neutrino 
radiation after their generation in the rotating plasma. Furthermore we assume that these 
particles can freely propagate along null geodesies even in the interior of the central object, 
and have only negligible energy losses in the interactions with the stellar matter. In the 
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following we determine the electron-positron energy deposition rate (EDR) produced due to 
the neutrino- anti neutrino pair annihilations into e~e + pairs along the rotation axis of the 
compact object, where the baryon contamination related to the matter content of the disk 
corona is minimal. 

For any point r = (r, 9, 0) of the u — v collision, the EDR per unit volume via electron- 
positron pair creation due to neutrino-anti neutrino annihilation is given by 

q(r) = = J J f v (p v , r)f D (p 9 , r){a vV \v v - v p \e u e p }^^d 3 p u d 3 p p , (14) 

where p v and p p are the 3-momenta in the momentum space, f u and f p are the neutrino 
and anti-neutrino number densities, and v u , e u , v D and e p are the 3- velocities an d the energy 
of the colliding neutrino-antineutrino pairs, respectively ( jGoodman et al. 19871 ). The cross- 
section of the collision is denoted by a uP , and can be calculated by using the formula 

a vP = KG 2 F {e u e D - c 2 p u ■ p v ) = -KG 2 F p v ■ p 9 

(special relativity is considered), where K = (1 ±4sm 2 9w + 8 sin 4 9w)/G^, with sin 2 9w = 

0.23, and G 2 F = 5.29 x 10~ 44 cm 2 MeV -2 , respectively. By applying the decompositions 

p u = e v £l v and d 3 p u = E 2 u dE v dQ, v , with the solid angle vector Q, v pointing in the direction of 

p V) and with the assumption that the neutrino source in the disk emits particles isotropically, 

the integral in Eq. f JT4|) can be separated into an energy integral and an angular part. 

After evaluating the energy integral for fermions and inserting the general relativistic 

version of the scalar product of the 4-momenta given by Eq. ( !T3|) into the cross-section 

equation for the colliding vv pair, for a stationary and axially symmetric geometry we obtain 

the EDR at the rotational axis (9 = 0) per unit volume per unit time (dtdV = y/—gd A x) in 

the form (jAsano fc Fukuyama 200 ll ) 

dE 
dtdV 

21 ^(8)^*11,(8r.)e(r). (15) 

Here T e //(3r 9 ) is the effective neutrino temperature evaluated at 3 xr s = 6GM/c 2 , and ( is 
the Riemann zeta function. The angular part 0(r) of the EDR is given by 



q{r 



2cKG 2 F e(r) J J f v f p (e v + e p )el4de v de t 



0(r) = 1 



T e//( 3r 9) \9t^-9ttgu 

[•0m /•27T p2n 

x d9 p sin 9 P d<p v d<f p T^{R 
Je m Jo Jo 

xTq(_R p )[1 — sin9 u sin^p cos^^ — <p p ) 



M 

d9 v sin 9 V 
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— cos cos^p] 2 , (16) 

with the neutrino temperature T observed at infinity, and the angular boundaries 9 m and 
6m of the collision angle for the (anti)neutrinos arriving form the disk edges at _R; n and R ou t, 
respectively. The expression in parenthesis is the redshift factor for massless particles. 

Integrating this expression over the azimuthal collision angles <p u and <£>p of the neutrino- 
antineutrino pairs yields the formula 
2tt 2 

— I ... .. . I " I 



4 r 

' — - I ii \ 

6(r) 



2 / dfi u T (fi v ) / c?/i p T (/ip) 

*f LLm. " Ltm. 



T eff(3r g ) V 9%, - 9tt94>4> e=0 , 
+ / dfi v (l - vl)T^ v ) / d^(l - /4)T*0 

+2 / d\i v \i y T Q (fj, u ) / dfi p fi P T (// p ) 



-4 



/Am " /Am 

MA/ _ /'MM 



5 /"^ M 4 

dfj Ju fj JU T (fi u ) / d(jLpfjbpT (fj,p) 



(17) 



where = cos^, /x p = cos^p, \x m = cos9 m and \im = cos 9m- 

In order to determine the boundary angles 6 m and #m (or, in turn, [i m and /xm), we 
traced the null geodesies back from a given point on the rotational axis to its intersection 
with the equatorial plane. In the ray tracing code we supplemented the coordinates of r 
and 9 with their momenta p T and pg, and used the dynamical equations for the extended 
set of variables. We present some details of this method in Appendix B. For a given set of 
canonical data of the final point at the rotational axis the code derived the canonical data 
of the initial point in the equatorial plane. That is, for any radial coordinate r along the 
axis, the code can provide the values of pg and p r for which null geodesies starting from the 
equatorial plane at a given radius R will hit the axis at r. The collision angle 9 V is determined 
by the scalar product of the normal vector pointing into the direction of the 3-momentum 
Pu = (iPv) r , (Pu) e , ij>uY) °f the neutrino passing through the axis, and the normal vector 
e z = {l/\fg~rr, 0, 0) of the rotational axis, 



cos 6^ = — — — = , (18) 

l P "l \fg rr [(Pv)r/grr + {Pufg/ Qee] 

where we have used the expression p u ■ e z = g r r(Pu) r (e z ) r = -Jg^ipuY , and the fact that the 
norm of the 3-momentum is reduced to \p u \ = \jg rr {Pv) 2 r + g ee (Pu)g, since (p„)<£ = L u is zero 
for the particles reaching the rotational axis. 

We can express the square of p r from the geodesic equation (j2J) as 
(Pv)r = gj4^o[cf (Pu) 2 g = gjiuJo/c) 2 /- 1 - ( Pu )j ^ 
gt<p ~ 9u9h> 9ee 9ee 
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Here the constant of motion loq for the particle energy is arbitrary (< c 2 ), provided the mass 
shell condition p v ■ p v — for null particles holds. Combining Eqs. ( lT8i) and ( fl9l) together 
we obtain a relation between the collision angle and pg: 



cos = ± 
or 

sin 2 9 V 



\ 



-i 2 



v 6 



{Pu)e 
u /c_ 



f(r) 



(20) 



(21) 

A/(r 2 + a 2 ), with 



u /c\ gl(r,0) ' 

From the latter formula in the Kerr spacetime we have sin9 u = p v 
the horizon function A = r 2 + a? — 2Mr, and p u = c(p u )g/uo = cJ/uq at the rotational axis 
( lAsano fc Fukuyama 20011 ) . We used the relation ( 1201) to calculate 9 V for any particle traced 
back from a given point at r on the rotational axis, with a given value of pg, to its starting 
point at R in the equatorial plane. For R = R in and R = R ou t, Eq- (1201 provides the values 
of 9 m , 9m (o, equivalently, /i m and /xm), for any value of r. 

For neutrinos propagating from the accretion disk with zero angular momentum in the 
direction of the disk rotation, the energy is redshifted by the gravitational potential, and 
Lorenz boosted due to the disk motion measured in the locally non-rotating frame (LNRF) 
of the spacetime. Then the observed neutrino temperature can be written as 

rim - \/ 9 ^~ 9tt9 ^ T (m 

I 0[R) - I eff[R) , 



(22) 



where T = (1 — v 2 /c 2 ) 1 / 2 , and the square of the disk velocity measured in the LNRF is 

— (tt K - uf 

9t4> - 9tt9<p<t> 

with the angular frequency Qk of the Keplerian disk, and the frame dragging frequency u 
of the spacetime given by 



-9t<i>,r + \J (gtcf>,r) 2 — 9tt, 



UJ 



dt g<f><f>,r gcjxj> 

respectively. If we solve Eq. (120]) for pg at a given value of 9 U , and use the result to prescribe 

the canonical data of any final point on the axis, then we can trace null particles back to 

their initial point in the equatorial plane. Hence we can determine the radial coordinate R 

measured along the disk at which To has to be evaluated in the integral given by Eq. ( TT7l) 

for any \i v = arccos 9 V . 

The energy deposition rate per unit volume for a distant observer is given by Eq. ffT5l) . 

with an integral over the volume dV = y/—gdrd9d<f). Since we estimate the EDR only along 
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the rotational axis over r min and r max with 9 — > 0, it is proportional to the integral 



dEr 



dt 



d 

2 d9 



2tt rir/2 rr max dEft 



dtdV 



-gdrd9d<p 



21vr 5 C(5)^ffc 9 T e 9 // (3r ff )r 2 



h 6 c 5 e//v 9 7 g 

where the dimensionless quantity G(r) is defined by 



G(r)dr 



G{r) 



-g(r,9) 



sin 9 



6(r) . 



(23) 



(24) 



We note that since g ~ sin 2 9 holds for the metric determinant of stationary and axially 
symmetric spacetimes in both the curvature and isotropic coordinate systems, dy/—g/d9\g =0 
can be written in the form {yj—gj sin9)(dsm9/d9)e = o = yf—gj sin#, which has a finite value 
at 9 = 0. For Kerr the metric the cofacto r of Q(r) in the Eg. ([24")) i s just (r 2 + a 2 )/r 2 , and 
G reduces to its definition introduced in (lAsano &: Fukuyama 20011 1 . 



4 EQUATIONS OF STATE AND STELLAR MODELS 

In order to obtain a consistent and realistic physical description of the rotating general 
relativistic neutron and quark stars, as a first step we have to adopt the equations of state 
for the dense neutron and quark matter, respectively. In the present study we consider the 
following equations of state for neutron and quark matter: 



1) Akmal-Pandharipande-Ravenhall (APR) EOS f lAkmal et al. 19981 ) . EOS APR has 



been obtained by using the variational chain summation methods and the Argonne vi$ two- 
nucleon interaction. Boost corrections to the two-nucleon interaction, which give the leading 
relativistic effect of order (v/c) 2 , as well as three-nucleon interactions, are also included in 
the nuclear Hamiltonian. The density range is from 2 x 10 14 g/cm 3 to 2.6 x 10 15 g/cm 3 . The 
maximum mass limit in the static case for this EOS is 2.20M^. We jo in this equation of 
state to the composit e BBP fe/c 2 > 4.3 x 10 n g/cm 3 ) (IBavm et al. 1971aT) - BPS (10 4 g/cm 3 



< 4.3 x 10 n g/cm 3 ) f lBavm et al. 1971bh - FMT (e/c 2 < 10 4 g/cm 3 ) f lFevnman et al. 19491 ) 



equations of state, respectively. 



2) Douchin-Haensel (DH) EOS ( IDouchin fc Haensel 20011 ). EOS DH is an equation of 



state of the neutron star matter, describing both the neutron star crust and the liquid core. 
It is based on the effective nuclear interaction SLy of the Skyrme type, which is particularly 
suitable for the application to the calculation of the properties of very neutron rich matter. 
The structure of the crust, and its EOS, is calculated in the zero temperature approximation, 
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and under the assumption of the ground state composition. The EOS of the liquid core is 
calculated assuming (minimal) npefi composition. The density range is from 3.49 x 10 11 
g/cm 3 to 4.04 x 10 15 g/cm 3 . The minimum and maximum masses of the static neutron stars 
for this EOS are O.O94M and 2.O5M , respectively. 

3) Shen- Toki-Oyamatsu-Sumiyoshi (STOS) EOS (IShen et al. 1998T ). The STOS equation 
of state of nuclear matter is obtained by using the relativistic mean field theory with nonlin- 
ear a and u terms in a wide density and temperature range, with various proton fractions. 
The EOS was specifically designed for the use of supernova simulation and for the neu- 
tron star calculations. The Thomas-Fermi approximation is used to describe inhomogeneous 
matter, where heavy nuclei are formed together with free nucleon gas. The temperature is 
mentioned for each STOS equation of state, so that, for example, STOS0 represents the 
STOS EOS for T = MeV. For the proton fraction we chose the value Y p = 10~ 2 in order 
to avoid the negative pressure regime for low baryon mass densities. 

4) Relativistic Mean Field (RMF) equations of state wi th isovector scalar mean fi eld cor- 
responding to the 5-meson- RMF soft and RMF stiff EOS (IKubis fc Kutschera 19971 ) . While 
the 5-meson mean field vanishes in symmetric nuclear matter, it can influence properties of 
asymmetric nuclear matter in neutron stars. The Relativistic mean field contribution due to 
the 5-field to the nuclear symmetry energy is negative. The energy per particle of neutron 
matter is then larger at high densities than the one with no 5-field included. Also, the proton 
fraction of /3- stable matter increases. Splitting of proton and neutron effective masses due 
to the 5-field can affect transport properties of neutron star matter. The equations of state 
can be parameterized by the coupling parameters C 2 = g 2 /m 2 , C 2 = g^/m^, b = bjg\ and 
c = c/g*, where m a and m w are the masses of the respective mesons, and b and c are the 
coefficients in the potential energy U (a) of the cr-field. The soft RMF EOS is parameterized 
by Cl = 1.582 fm 2 , C 2 = 1.019 fm 2 , b = -0.7188 and c = 6.563, while the stiff RMF EOS 
is parameterized by C 2 = 11.25 fm 2 , C 2 = 6.483 fm 2 , b = 0.003825 and c = 3.5 x HT 6 , 
respectively. 

5) Baldo-Bombaci-Burgio (BBB) EOS (IBaldo et al. 1997T ). The BBB EOS is an EOS 
for asymmetric nuclear matter, derived from the Brueckner-Bethe-Goldstone many-body 
theory with explicit three-body forces. Two EOS's are obtained, one corresponding to the 
Argonne AV14 (BBBAV14), and the other to the Paris two-body nuclear force (BBBParis), 
implemented by the Urbana model for the three-body force. The maximum static mass con- 
figurations are M max = 1.8M Q and M max = 1.94M Q when the AV14 and Paris interactions 
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are used, respectively. The onset of direct Urea processes occurs at densities n ^ 0.65 fm -3 
for the AV14 potential and n ^ 0.54 fm -3 for the Paris potential. The comparison with 
other microscopic models for the EOS shows noticeable differences. The density range is 
from 1.35 x 10 14 g/cm 3 to 3.507 x 10 15 g/cm 3 . 



6) Bag model equation of s tate for quark matter (Q) EOS (lltoh 1970 



Bodmer 1971 



Witten 1984 



Cheng et al. 19981 ). For the description of the quark matter we adopt first a 



simple phenomenological description, based on the MIT bag model equation of state, in 
which the pressure p is related to the energy density p by 



p=-(p-AB) c 2 



(25) 



where B is the difference between the energy density of the perturbative and non-perturbative 
QCD vacuum (the bag constant), with the value 4B = 4.2 x 10 14 g/cm 3 . 

7) It is generally agreed today that the color-flavor-locked (CFL) state is likely to be the 
ground s tate of matter, at l east for asymptotic densities, and even if the quark masses ar e 
unequal f Alford et al. 1999 : Rapp et al. 2000 : Horvath fc Lugones 20pi Alford et al. 2007 ). 



Moreover, the equal number of flavors is enforced by symmetry, and electrons are absent, 
since the mixture is automatically neutral. By assuming that the mass m s of the s quark 
is not large as compared to the chemical potential p, the thermodynamical pot ential of the 
quark matter in CFL phase can be approximated as ( ILugones Horvath 20021 1 



n 



CFL 



3/i 4 3m 



4vr 2 



+ 



4vr 2 



121n(m a /2/x) 
32vr 2 



A V + B, 



7T 



(26) 



where A is the gap energy. With the use of this expression the pressure P of the quark 
matter in the CFL phase can be ob tained as an explicit function of the energy density e in 
the form (ILugones fe Horvath 2002l ) 

m 2 5 2 



1 , x 2A 2 5 2 

P= -(e-4B +— 5- 

3 7l Z 

where 



2tt 2 



4 



(27) 



(28) 



<T = -a + ya 2 + -7r 2 (e-B), 

and a = — m 2 /6 + 2A 2 /3. In the following the value of the gap energy A considered in each 
case will be also mentioned for the CFL equation of state, so that, for example, CFL200 
represents the CFL EOS with A = 200. For the bag constant B we adopt the value AB = 
4.2 x 10 14 g/cm 3 , while for the mass of the strange quark we take the value m s = 150 MeV. 
The pressure-density relation is presented for the considered equations of state in Fig. [TJ 
In order to calculate the equilibrium configurations of the rotating neutron and quark 
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Figure 1. Pressure as a function of density (in a logarithmic scale) for the equations of state DH, RMF soft, RMF stiff, STOS 
0. STOS 0.5, STOS 1, BBBAV14, BBBParis, APR, Q, CFL150, and CFL300, respectively. 

stars with the EOS's presented here we use the RNS code, as introduced in Stergioulas & 
Friedman (1995), and discussed in details in Sterg ioulas (2003). This c ode was used for the 
study of different models of rotatin g neutron stars (jNozawa et al. 19981 ). and for the study of 
the rapidly rotating strange stars (jStergioulas et al. 19991 ). The metric outside the rotating 
compact general relativistic stars can be described, in quasi-isotropic coordinates, as 

- udt) 2 , 



ds 2 = -e~< +p dt 2 + e 2a (df 2 + f 2 d0 2 ) + e 7 "^ 2 ~'- 2 



r sin 



(29) 



where the metric potentials 7, p, a and the angular velocity of the stellar fluid relative to 
the local inertial frame u are all functions of the quasi-isotropic radial coordinate f and of 
the poloidal angle 9. The RNS code computes the metric functions in a quasi-spheroidal 
coordinate system and write them as functions of the compactified dimensionless distance 
s = fj (f + f e ), where f e is the equatorial radius of the star, and the cosine of the poloidal 
angle. The quasi-spheroidal radial coordinate is converted into curvature coordinate r ac- 
cording to the equation r = f exp {[7(r, 9) — p(f, 9)} /2}. 



5 THE ISOTHERMAL THIN DISK 

With the assumption that the surface of the thin accretion disk is in an isothermal state, 
(which is definitely a crude approximation), the effective temperature T e jf of the d isk surface 



can be set to a constant over the entire disk. Following ( lAsano fc Fukuyama 20011 ) . we use a 
constant effective temperature measured at 3r g , T e ff = T e ff(3r g ), with T e ff in the range of 10 
to 20 MeV. In this case the radial dependence of the neutrino temperature T is determined 
only by the gravitational redshift, and the Lorentz boost, as can be seen in Eq. (122]) . Since 
the neutrinosphere either does not exist (if the central body is not compact enough), or its 
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equatorial radius is smaller than the stellar radius, for the configurations of the neutron and 
quark star models pr esented here, the en tire surface of the accretion disk can be considered 
as a neutrino source (jKovacs et al. 20101 ). 

As it has been showed before, if black holes spin up the energy de position via the neutrino 



pair annihilation i s increasing at the rotational axis of the black hole (jAsano fc Fukuyama 2001 



Miller et al. 20031 ) . This phenomenon is a consequence of the decreasing radius of the marginally 



stable orbit, i.e., of the increasing surface area of the accretion disk with the increasing spin. 
However, this enhancement is moderated by the redshift of the neutrinos emerging from the 
innermost region of the disk, which causes a smaller and smaller increment in the EDR as 
the spin parameter is approaching one. The Doppler shift due to the disk rotation has also to 
be taken into account when we calculate the net temperature effect on the axial distribution 
of the EDR. 

If the central body is a neutron or quark star, then the situation is even more complicated. 
In many cases there is no marginally stable orbit outside the star and the inner edge of the 
accretion disk is located at the stellar surface. Then the proportionality between the area of 
the disk surface and the spin is no longer guaranteed, and it indeed breaks down at higher 
rotational frequencies for some physical configurations of stars, and for some EOS types 
where inner disk edge touches the stellar surface. 

In Table [1] we present the relation between the spin parameter, the radius of the inner 
edge of the disk, and the enhancement in the EDR. Here we have integrated the EDR given 
by Eq. (1231 over r = 3 — 10r 9 for each configuration, EDR = J^ s (dE /dt)g =0 dr, and we 
have computed AEDR=(EDR-EDR S )/EDR S , where EDR S is derived from the Schwarzschild 
geometry. The cutoff radius R ou t of the accretion disk was set to 10r 9 . In the calculations 
we used the metric provided by the RNS code for the total mass M = 1.8M , and different 
angular velocities between 5 x 10 2 s _1 and 5 x 10 3 s _1 . The tables containing the physical 
parameters of these stellar configurations can be found in Appendix C. 

For the first four stellar models presented in Table HJ namely those with DH, APR, 
BBBAV14 and BBBParis type EOS's, the EDR is increasing as the disk surface is increasing 
(Ri n is decreasing) with the increasing spin of the star. Since at the rotation frequency of 
5 x 10 3 s _1 the value of the spin parameter is only about 0.3 — 0.35, the increase in the EDR 
as compared to the EDR obtained for a static star or black hole is only ~ 7%. For a black 
hole spinning up to a* = 0.99, this enhancement is already 30% (jAsano fc Fukuyama 20011 ) 

We also present the axial distribution of the EDR for these stellar models in Fig. [2j 
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Figure 2. The axial profile of the EDR characterized by the dimensionless quantity G(r) as a function of f! for an isothermal 
accretion disk and neutron star models with DH (top left hand side), APR (top right hand side), BBBAV14 (bottom left hand 
side) and BBBParis (bottom right hand side) type EOS. The total stellar mass is 1.8Mq. For comparison we also plotted G(r) 
for a static spacetime geometry (black solid curve). 



where we have plotted the quantity G(r), defined in Eq. (124|) . as a function of the curva- 
ture coordinate r. The graphs show that the proportionality between the maximum of the 
deposition rate and the spin parameter of the star is also preserved: the faster the neutron 
star rotates, the higher the maximal value of the EDR along the axis. Nevertheless, these 
maxima shift closer to the stellar surface for higher angular velocities, which involves the 
baryon contamination problem. Besides the uncertainties arisen from baryon contamination, 
one must be careful with the interpretation of these results, since the application of a for- 
malism developed for black holes may produce considerable errors at higher values of Q. 
As the error due to the separation of the geodesic equation is increasing with the higher 
rotational speed of the central object, the calculations produce less reliable results for high 
Q. Hence it is not clear that the moderate increase in the EDR at higher angular frequencies 
is the result of the effect of the redshift of the neutrino/anti neutrino energy alone, or it is 
caused by the underestimation of the real values in the approximation. The latter must at 



18 Kovdcs, Cheng and Harko 



EOS 


n [10 3 s- 1 ] 


0.5 


1.0 


1.5 


2.0 


2.5 


3.0 


3.5 


4.0 


4.5 


5.0 


DH 


CI* 


0.03 


0.07 


0.10 


0.13 


0.17 


0.21 


0.24 


0.27 


0.31 


0.35 




■^^in 1 ^ Q 


2.94 


2.89 


2.85 


2.79 


2.75 


2.70 


2.67 


2.63 


2.60 


2.56 




AEDR (isothermal) [%] 


1.11 


1.95 


2.59 


3.38 


4.17 


5.20 


5.98 


6.19 


6.96 


7.56 




AEDR (T e ff(R)) [%] 


12.8 


25.0 


37.0 


51.9 


66.9 


85.9 


100.6 


115.8 


133.3 


149.6 


APR 


CI* 






0.10 


0.14 


0.17 


0.20 


0.24 


0.27 


0.31 


0.35 










2.84 


2.79 


2.75 


2.71 


2.66 


2.63 


2.59 


2.56 




AEDR (isothermal) [%] 


- 


- 


1.51 


3.35 


4.16 


5.00 


5.82 


6.44 


6.70 


7.59 




AEDR (T e ff(R)) [%] 


- 


- 


36.5 


53.1 


66.6 


83.3 


102.0 


115.7 


134.6 


149.1 


BBBAV14 




0.03 


0.05 


0.08 


0.11 


0.13 


0.16 


0.19 


0.23 


0.26 


0.29 






2.95 


2.91 


2.88 


2.83 


2.79 


2.75 


2.71 


2.66 


2.62 


2.58 




AEDR (isothermal) [%] 


0.19 


1.12 


1.81 


2.73 


3.54 


4.37 


4.91 


5.62 


6.30 


7.14 




AEDR (T e ff(R)) [%] 


9.8 


19.4 


28.5 


41.2 


53.5 


68.8 


83.7 


102.2 


120.2 


143.2 


BBBParis 


(I* 




0.06 


0.09 


0.12 


0.15 


0.18 


0.21 


0.24 


0.27 


0.31 








2.90 


2.86 


2.81 


2.77 


2.73 


2.69 


2.64 


2.61 


2.57 




AEDR (isothermal) [%] 


- 


1.71 


2.49 


3.38 


4.11 


4.84 


5.57 


6.39 


6.47 


6.94 




AEDR (T e ff(R)) [%] 


- 


21.0 


32.8 


47.3 


60.2 


74.4 


90.0 


111.5 


124.5 


144.7 


RMF stiff 


a* 


0.05 


0.09 


0.14 


0.20 


0.25 


0.30 


0.36 


0.43 


0.51 


0.63 




^^XTl / ^* Q 


2.92 


2.87 


2.81 


2.77 


2.75 


2.63 


2.67 


2.75 


2.83 


2.96 




AEDR (isothermal) [%] 


1.65 


2.62 


3.65 


4.37 


4.88 


6.03 


5.72 


4.50 


3.29 


2.09 




AEDR (T e ff(R)) [%] 


16.8 


32.0 


48.0 


62.4 


70.0 


114.6 


96.8 


68.7 


40.8 


6.5 


STOSO 


a* 


0.06 


0.12 


0.17 


0.24 


0.30 


0.38 


0.46 


0.56 


0.69 


0.86 




-^^171 f 9 


2.92 


2.87 


2.88 


2.90 


2.96 


3.00 


3.10 


3.21 


3.43 


3.84 




AEDR (isothermal) [%] 


1.88 


2.69 


3.05 


2.96 


1.96 


1.65 


-0.41 


-1.49 


-6.08 


-12.36 




AEDR (T e ff(R)) [%] 


21.7 


36.2 


34.7 


27.3 


14.8 


4.5 


-15.0 


-31.2 


- 56.6 


-80.0 


Q 


0* 






0.11 


0.15 


0.19 


0.23 


0.28 


0.32 


0.37 


0.41 










2.83 


2.78 


2.74 


2.71 


2.67 


2.67 


2.67 


2.67 




AEDR (isothermal) [%] 






2.99 


4.01 


4.54 


5.11 


5.48 


5.82 


5.92 


6.37 




AEDR (T e ff(R)) [%] 






40.2 


55.8 


70.5 


84.1 


96.3 


101.6 


102.5 


100.2 


CFL150 


a* 




0.08 


0.13 


0.17 


0.21 


0.26 


0.31 


0.36 


0.41 


0.47 








2.88 


2.82 


2.77 


2.73 


2.71 


2.69 


2.70 


2.73 


2.79 




AEDR (isothermal) [%] 




2.15 


3.07 


4.10 


4.59 


5.06 


5.30 


5.47 


5.20 


5.06 




AEDR (T e ff(R)) [%] 




28.7 


44.1 


60.3 


72.5 


82.3 


88.6 


88.3 


78.0 


59.0 



Table 1. The spin parameter a* of the star, the inner edge Ri„ of the accretion disk, and the variation in the EDR integrated 
over the rotational axis for neutron and quark star models, with a total mass of 1.8Mq and angular velocities from 5 X 10 2 s~ 1 
to 5 X 10 3 s -1 . Here AEDR = (EDR - EDR S )/EDR S and EDR S is derived for the static case (fi = 0). The variation in EDR 
is calculated for both isothermal disks (T e yy=const) and disks with the temperature profile T e ff(R) = T e ff X 3r g /R. 



least not dominate over the real physical effects, as far as we consider configurations with 
lower values of the spin parameter (see Appendix A). 

From the off-axis calculations of the spatial distribution of the EDR in the vicinity 
of Kerr black holes we already know that the increase in the rotational speed of the black 
hole enhances the e nergy deposition along both the equatorial plane, and the axis of rotation 
(jMiller et al. 20031 ). This may indicate that the behavior presented in Fig. |2]is in accordance 
with the previous resul ts on the EDR dist ribution along the equatorial plane for the same 
family of stellar EOSs ( IKovacs et al. 2010l ). It was shown that for the APR type EOS the 
maximum of the energy deposition from vv pair annihilation is increasing along the surface 
of the thin accretion disk as the star spins up. The other types of EOS shown in Fig. [2] 
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Figure 3. The axial profile of EDR characterized by the dimensionless quantity G(r) as a function of f2 for an isothermal disk 
and neutron star models with RMF stiff (top left hand side) and STOSO (top right hand side) type EOSs, and for quark stars 
with Q (bottom left hand side) and CFL 150 (bottom right hand side) type EOSs. 



exhibit similar behaviors. Therefore we can conclude that they follow the pattern found in 
the case of Kerr black holes. 

The neutron star models with RMF stiff and STOSO type EOS, together with the Q and 
CFL 150 type quark stars do not produce a simple relation between the spin and the EDR, 
as shown in the bottom four lines of Table [U In Fig. [3] we also present the dependence of 
the EDR on Q for these models. In the case of the neutron stars with RMF stiff and STOSO 
type EOSs there is some enhancement in the EDR for rotational velocities increasing up 
to ~ 3 x 10 3 s" 1 (RMF stiff) and ~ 1.5 x 10 3 s" 1 (STOSO), but at higher values of Q the 
maxima of the ERD starts to decrease. As the star rotates faster its surface reaches higher 
equatorial orbits, and reduces the surface area of the accretion disk by pushing the inner 
edge of the disk to higher radii. At rotational frequencies approaching 5 x 10 3 s -1 the amount 
of the energy deposition for the stellar model with RMF stiff type EOS reduces to the one 
derived in the static case. Since the stellar radius is already much greater than r ms = 3r g of 
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Figure 4. The axial profile of the EDR, characterized by the dimensionless quantities G(r) for an isothermal disk, and different 
neutron and quark star models with the total mass M = 1.8Mq and rotational velocities f2 = 2 X 10 3 s —1 (top left hand side 
panel), 3 X 10 3 s" 1 (top right hand side panel), 4 X 10 3 s _1 (bottom left hand side panel) and 5 X 10 3 s — 1 (bottom right hand 
side panel). 



the static case, we obtain even smaller values of the EDR for the STOS0 model then those 
derived for a static star (see the negative values in Tabled]). 

The quark star models with Q and CFL150 type EOS do not deviate so drastically from 
the spin-EDR relation characterizing the first family of EOSs and the black holes. In Fig. [3] 
we can see that the maximum of the EDR has an increase, turns over, and then decreases 
with increasing spin parameter. However, its reduction after the turn over is minimal, and 
it might be a result of the combination of several physical effects, such as the variation of 
the disk surface area, the redshift and the Lorenz boost, together with the artifacts caused 
by the break-down of our approximation for higher values of the spin parameter. The total 
EDR integrated along the rotational axis remains increasing for the Q type EOS, but the 
differences are small, and may be dominated by the uncertainties for fast stellar rotation. 
These results are still in good agreement with the fact that the amount of energy liberated in 
the vv pair annihilation along the equatorial plane is proportional to the angular velocity of 
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the quark star with Q type EOS, as far as we compar e the configurations for = 3x 10 3 s _1 
and 5 x 10 3 s _1 with each other (IKovacs et al. 20101 ) . 



In Fig. H] we also compare the EDR distributions calculated along the rotational axis 
for different EOS types at constant angular frequencies ranging from Q = 2 x 10 3 s _1 to 
5 x 10 3 s _1 , respectively. The plots show that the quark stars and the neutron stars with 
RMF stiff type EOS produce the highest values of the EDR at low rotational frequencies, 
although the differences are not considerable. This relation is similar to the one derived in 
the equatorial plane, where the EDR is higher for quark stars than for neutron stars, apart 
from those with RMF stiff and STOS0 EOS type (IKovacs et al. 20101 ) . The only difference 
is the behavior of the STOS0 model, which has somewhat lower values in the EDR along 
the rotational axis than other neutron star models d o. However, the ED R values in the 
equatorial plane were derived at = 5 x 10 3 s _1 in ( Kovacs et al. 201ol ). and we should 
compare them with the plots displaying configurations with high rotational velocities. 



The comparison shows that quark stars and neutron stars with RMF stiff and STOS0 
type EOS generate less energy along the rotational axis from the neutrino-anti neutrino 
annihilation than the neutron stars with other EOS types do at Q = 5 x 10 3 s _1 , as op- 
posed to the relations between the values of EDR calculated in the equatorial plane. In 
fact, there is no discrepancy between the two results, as the EDR at a given point R in 
the equat orial plane is calculated by integrating the angular part of Eq. (HM between R 
and R ou t (jSalmonson &: Wilson 19991 ). whereas the calculation of the EDR at a point r on 
the rotational axis always involves the integral of Eq. (Tbil) over the entire accretion disk 
(lAsano fc Fukuyama 20011 ) . As a result, the first method provides q for any value R > Ri n , 
which does not depend on R in (the latter one is only a lower boundary of the domain of q), 
whereas the EDR derived by the second method always depends on R in , i.e, the area of the 
disk surface. 



Nevertheless, the results presented here are already uncertain in the regime of high 
angular frequencies, and the EDR along the rotational axis should already be calculated for 
rapidly rotating stars by applying the fully 3D techniques of off-axis calculations. 
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Figure 5. The axial profile of the EDR characterized by the dimensionless quantity G(r) for T e ff oc R~ 1 and different angular 
velocities of the neutron star models with DH (top left hand side), APR (top right hand side), BBBAV14 (bottom left hand 
side) and BBBParis (bottom right hand side) type EOSs. (The total mass of the central object is I.SMq.) For comparison we 
have also plotted G(r) for a static spacetime geometry (black solid curve). 

6 THIN DISK WITH RADIATION IN THERMODYNAMIC AL 
EQUILIBRIUM 

In this Section we use for the calculation of the EDR along the rotational axis for various ro- 
tating compact objects a more realistic model of the thermal distribution of the neutrino flux 
emerging from the thin accretion disk. This model is based on the thermodynamical equilib- 
rium approximation of the state of the accretion disk, where the temperature gradient is given 
by a profile inversely proportional to the radius. Following Asano & Fukuyama (2001), we 
replace the constant temperature profile T e ff = T e jf(3r g ) with T e ff(R) = T e ff(3r g ) x 3r g /R, 
and we recalculate the EDR along the rotational axis for all the constant disk temperature 
configurations considered before. 

In Figs. [5] and [6] we present the axial distribution of the EDR derived for the accretion disk 
with temperature gradient. If we compare the results with those obtained in the isothermal 
case for the same T e ff(3r g ), we see that the EDR is reduced by a factor of 100 at low 
rotational frequencies of the central body. Even for the highest values of Q, the EDR in the 
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Figure 6. The axial profile of the EDR characterized by the dimensionless quantity G(r) for T e ff oc R~ 1 , and different angular 
velocities of the neutron star models with RMF stiff (top left hand side) and STOSO (top right hand side) type EOSs, and for 
quark stars with Q (bottom left hand side) and CFL 150 (bottom right hand side) type EOSs. 



isothermal case is about 25 times higher than the one calculated with temperature gradient. 
Only the dynamical range of the EDR became wider for the first group of EOS in Fig. |5j 
in the isothermal case its values at r ~ 3r g have a ~ 10% increase for Q = 5 x 10 3 s -1 , 
as compared with those for Q = 5 x 10 2 s _1 . The difference increases to a factor of 3 for 
the disk with the temperature gradient. Table [T] also contains the differences of the total 
EDR integrated over r = 3r g — 10r g for these configurations for a disk with the /^-dependent 
temperature profile. The differences in the total EDR compared with the static case are 
from ~ 10% to ~ 150%, which span a much wider range than those for the isothermal disk 
do, between ~ 1 and ~ 7%. Although in the models with a disk temperature gradient the 
efficiency of the energy deposition compared with the static case increases in a much higher 
rate with increasing Q, the comparison of Figs. [2] and [5] shows that even at high rotational 
velocities the EDR is much smaller for T e ff oc than the one calculated for isothermal 
disks. The plots also demonstrate that the EDR decreases more rapidly with r for the latter 
model than it does for isothermal disks, especially at higher angular velocities. This means 
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that the axial distribution of the EDR is more concentrated to the region close to the star, 
and the energy is produced here mainly by the neutrinos and anti - neutrinos emitted from 
the innermost region of the disk. For the temperature profile T e ff oc Rr 1 , this region is 
hot, but the neutrinos leaving the rapidly rotating disk surface suffer a large Doppler shift, 
which, together with the large gravitational redshift, prevents them from depositing much 
energy at higher radii. 

For the second group of the EOS types shown in Fig. [HI the EDR still preserves its non- 
linear dependence on the rotational velocity of the central body. Over fl = 3 x 10 3 s -1 for 
the RMF stiff model and over 10 3 s _1 for the STOSO EOS type their relation is inverted, 
and the EDR is approaching to, or even takes smaller values, than the one obtained in the 
static configuration. This trend is also present in the differences of the total EDR integrated 
along the axis, as can be seen in Table [TJ where AEDR has a considerable reduction at 
higher values of Q, as compared with its maxima. 

The quark stars follow the similar trend, but with moderate deviations, as in the isother- 
mal case. For higher rotational speeds in Table Q] the proportionality between the spin and 
the EDR is inverted, and we obtain a slight decrease in the total EDR as compared with 
the cases of slower rotation at 3.5 — 4 x 10 3 s -1 . 

In Fig. [7] we present G(r) for accretion disks with a T e ff oc R~~ l temperature profiles, 
rotating around stars with different EOSs types, but with fixed values fl These plots also 
show the same trends as those found in the case of the isothermal disk: at lower rotational 
frequencies the EDR is somewhat higher for the quark stars than the one derived for neutron 
stars - apart from the stellar model with RMF stiff type EOS - but with increasing Q its 
value for the Q and the CFL150 type EOS drops below those of the stellar models with DH, 
APR, BBBAV14, and BBBParis EOSs, respectively. The only neutron stars producing less 
EDR at higher frequencies than the one obtained for the quark stars are the models with 
RMF stiff and STOSO type EOS. 

7 DISCUSSIONS AND FINAL REMARKS 

In the present paper we have considered the energy deposition rate from the neutrino- 
antineutrino annihilation process on the rotation axis of rotating neutron and quark stars. 
The energy deposition rate has been obtained numerically for several equations of state of 
the neutron matter, and for two types of quark stars, respectively. All the general relativistic 
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Figure 7. The axial profiles of the EDR characterized by the dimensionless quantities G(r) for an isothermal disk, and different 
neutron and quark star models, with the total mass of the central object M = 1.8A/q, and rotational velocities f2 = 2 X 10 3 s _1 
(top left hand side panel), 3 X 10 3 s _1 (top right hand side panel), 4 X 10 3 s _1 (bottom left hand side panel) and 5 X 10 3 s _1 
(bottom right hand side panel). 



correction factors, related to this process, can be obtained from the metric of the central 
compact object. In the present paper we have neglected the possible effects of the stellar 
magnetic field on the external geometry of the spacetime around the compact general rel- 
ativistic effects, as well as on the motion of the neutrinos. Generally, we assume that the 
magnetic energy of the disk, B 2 /8ir is much smaller than the the kinetic energy 

density of the disk, B 2 /8tt << pV^/2, where Vk is the Keplerian velocity. Hence we neglect 
the possible influence of the magnetic field on the temperature profile of the disk. Due to the 
differences in the space-time structure, the quark stars present some important differences 
with respect to the energy deposition rate, as compared to the neutron stars. As a general 
result we have found that there is a general correlation between the energy deposition rate 
and the spin parameter of the star. There is also a strong dependence between the temper- 
ature profile of the disk and the EDR, and at high rotational velocities the EDR is much 
smaller for the model with T e ff oc R^ 1 than the one calculated for the constant temperature 
(isothermal) disks. Our study is based on the main assumption of the separability of the 
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geodesic equations of motion for the neutrinos. Even that from a purely theoretical point of 
view the geodesic and the Hamilton- Jacobi equations cannot be separated, from a practical- 
computational point of view we have found that such a procedure can work, at least for low 
angular speeds, with a good precision, with the error introduced by this approximation not 
exceeding 1%. 

In order to compare the properties of the quark and neutron stars we have chosen a 
common mass value M = 1.8M . Generally, quark stars may have masses greater than 
1.8M Q . However, for some equations of state of the neutron matter the maximum allowable 
mass could be smaller than the maximum mass of quark stars. Therefore we have chosen 
this particular value of the stellar mass in order to be able to compare sequences of quark 
stars and neutron stars for a large number of equations of state. 

A possible astrophysical application of our results could be in the explanation of the 
physical processes that lead to the formation of the Gamma- Ray Bursts (GRBs). The so- 
called fireball model can basically explain the observational facts well, and thus it is strongly 
favored, and widely accepted today (for recent reviews on GRBs see Meszaros (2006) and 
Zhang (2007), respectively). Newborn neutron stars, surrounded by hyper-accreting and 
neutrino-cooled disks may exist in some gamma-ray bursts and/or supernovae. As compared 
with the black hole disk, the neutrino annihilation luminosity above the neutron star disk 
is higher (IZhang fc Dai 20091 ). Although a heavily mass-loaded outflow from the neutron 
star surface at early times of neutron star formation prevents the outflow material from 
being accelerated to a high bulk Lorentz factor, if the disk accretion rate and the neutrino 
emission luminosity from the surface boundary layer are sufficiently high, an energetic ultra- 
relativistic jet via neutrino annihilation can be produced above the stellar polar region at 
late times. The energy deposition along the rotation axis will form a narrow cone jet, which 
encounters much less material, and hence it can maintain a large Lorentz factor. 

The neutrino emission rate is also strongly dependent on the temperature. In the stan- 
dard models of GRBs it is assumed that the central object is surrounded by a degenerate 
accretion di sk, which allo ws super-Eddington accretion rates, of the order of one solar mass 
per second (IZhang 20071 ). If the central compact object is a neutron or a quark star, such 
super-Eddington accretion rates can maintain the compact object at very high tempera- 
tures, and thus allowing very high neutrino luminosities, as well as a high rate of electron- 
positron pair production. The neutrino temperature can be estimated by assuming that 
the accretion power r]Mc 2 , where r\ is the efficiency of the energy conversion and M is 
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the accretion rate, is equal to the radiation power 4irR 2 aT^, which gives the temperature 
as T v = [f]Mc 2 1 '4:7iR 2 a^ ^ . By taking n = 0.1, R — 10 6 cm, and an accretion rate of 
M = 1Mq/10 s, we obtain T v = 7.14 x 10 10 K, a temperature which is of the order of MeV. 
Therefore, if the accretion disk is fed at a high rate, like, for example, by the fallback material 
after a supernova explosion, a high neutrino-antineutrino emission rate can be maintained, 
and this could explain some of the basic properties of the GRB phenomenon. 

The neutrino annihilation and the electron-positron pair production also plays an essen- 
tial role in the astrophysical processes associated with the phase transitions that could take 
place inside neutron stars. For example, the sudden phase transition from normal nuclear 
matter to a mixed phase of quark and nuclear matter induces temperature and density os- 
cillations at the neutrinosp here. Consequently , pulsating neutrino/antineutrino and leptonic 
pair fluxes will be emitted ( ICheng et al. 20091 ). During this stage several mass ejecta can be 
ejected from the stellar surface by the neutrinos and antineutrinos. These ejecta can be fur- 
ther accelerated to relativistic speeds by the electron/positron pairs, created by the neutrino 
and antineutrino annihilation outside the stellar surface. In order to produce the Gamma- 
Ray Bursts, a high neutrino emission rate is necessary. On the other hand, it is important to 
note that electron-positron pairs can deposit energy much more efficiently than the neutri- 
nos, and the dominant ene rgy deposition process is the neutrino-antineutrino annihilation 
process (ICheng et al. 20091 ). In fact most pairs are created outside the star. A large fraction 
of the neutrino energy, will be absorbed by the matter very near the stellar surface. When 
this amount of energy exceeds the gravitational binding energy, some mass near the stellar 
surface will be ejected, and this mass will be further accelerated by absorbing pairs created 
from the neutrino and antineutrino annihilation processe s outside the star. This process may 
be a possible mechanism for short Gamma- Ray Bursts (ICheng et al. 20091 ). 
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Figure Al. The Ag_|_/{g_|_) vs. s relation for neutron star models with various EOS types, such as DH, shown in Table ICTl 
(top left hand side plot), APR , shown in Table [C2] (top right hand side plot), BBBAV14, shown in Table IC3l (bottom left hand 
side plot) and BBBParis shown in Table IC4I (bottom right hand side plot). 



APPENDIX A: THE SEPARABILITY OF THE GEODESIC EQUATION 

The condition for the separability of the geodesic equation Eq. <^ for L = is that the 
functions / and g + , given by Eqs. (J3J) and (JSJ), respectively, must depend only on the radial 
coordinate r. This condition does not hold for general stationary and axially symmetric 
spacetimes, since the metric components are also functions of the polar angle. Now we ex- 
amine if we can use the approximation of neglecting the angular dependence of the functions 
/ and g + for the spacetime geometry of the neutron and quark stars which are considered in 
the present paper. If we insert the components g a/ 3 of the metric fl29|) . derived by the RNS 
code for these types of stars into Eqs. (J3J) and jSJ), we obtain f{f) = f 2 , and an expression 
for g + which still depends on both f and 6. Since / is the square of the radial coordinate 
only, the question is the angular dependence of g + . We consider the quantity Ag + / (g + ) as a 
function of the dimensionless radial coordinate s = f/(f + f e ) (compactified in the domain 
[0,1]) and angle variable \x = cos 8. Here Ag + is the maximal difference in the values of 
g + , and (g + ) is the average of g + for a given s over the entire domain [0,1] of \i. If g + has 
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Figure A2. The Ap+/(g+) vs. s relation for neutron and quark star models with various EOS types, such as RMF stiff, 
shown in Table IC5l (top left hand side plot), STOS0 , shown in Table ICol ftop right hand side plot), Q, shown in Table IC7l 
(bottom left hand side plot) and CFL150 shown in Table IC8I (bottom right hand side plot). 



no angular dependence, as in the case of the Kerr spacetime (in both Boyer-Lyndquist and 
quasi-isotropic coordinates), this quantity is just zero. Its deviation from zero characterizes 
the limits of the validity of our approximation. In Figs. IA1I and IA2I we plot this function 
versus the coordinate s for all the stellar models we have studied in this work. 

By definition of the coordinate s, the stellar surface in the equatorial plane is located 
at s — 0.5 and still remains close to this value at higher latitude. Therefore, the half of the 
right hand side of each plot represents the region of the spacetime outside the star. The 
accretion disk is rotating in this region, and the null geodesies emanating from the disk 
and intersecting the rotational axis at radii higher than 3r g mostly remain located in this 
region as well. The trajectories originating at the inner disk edge might run inside the star, 
but we consider only those particles which hit the axis at 3r g or higher radii. However, this 
minimal value is high enough to keep these trajectories close to the stellar surface (s < 0.5). 
In the case of the low mass stellar models with M = 1.8M , Ag + /(g + ) is under 0.5% for 
Q = 2 x 10 3 s _1 , while it is more the 1% for angular velocities higher than 4 x 10 3 s _1 . Then 
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on the one hand we do not expect that the separation of Eq. (j2J) causes large errors in the 
computation of scalar product p ■ p for slowly rotating objects, with Q ~ 2 x 10 3 s _1 . On 
the other hand, the results derived for Q > 4 x 10 3 s _1 have an increasing uncertainty as the 
error propagation grows due to the non-separability of Eq. (j2J), containing g + (r,9), which 
does depend on 9. 



APPENDIX B: THE RAY TRACING ALGORITHM 

The Hamilton- Jacobi equation (|2J) can be written as 

^r 2 + e 2 = geeV(u ,L), (Bl) 

9rr 

where the effective potential is given by 

y, L ^ _ 9<m>(uq/c) 2 + 2g t4> (ui /c)L + g tt L 2 

9t<f> ~ 9tt9H> 

This is generally not a separable differential equation for the coordinates r and 9. For given 
initial values of the coordinates at some time t = Tq, and for fixed constants of motion uq and 
L, we can, at least numerically, integrate this equation, to obtain the world-line projected 
onto the r — 9 plane of a particle moving in the effective potential V. The only technical 
difficulty in the integration may stem from the non-separability of Eq. (IBlj) . However, this 
can be easily circumvented by introducing the momenta p r and pp a nd by replacing Eq. (IBlj) 



with their canonical equations of motion (IPerez-Giz fc Levin 20091 ). Then the complete set 
of the evolution equations is 

f = (B2) 

9rr 

e = —, (B3) 



Pr = UtYt + tfri-V,], (B4) 

Z V 9r 



1 ( 9rr,e 2 , 9ee,e 2 t/ 1 

- 1 ~^-Pr + ~^-Pe -V,e), (B5) 



2 V 9 2 r 

together with the decoupled geodesic equation for the coordinates t and 0. This is a system 
of ordinary differential equations, with the unknown functions r(r), 9{t) p r (r) and pe(T~), 
respectively. The parameters appearing in the system are the constants of motion p t = —uq/c 
and p,!) = L. If the initial values of the unknown functions and of the parameters are given, 
so that the mass-shell condition p ■ p = holds (a constraint imposed on the initial data), 
then the system (1B2j) - (1B5j) can be integrated numerically. We have applied the fourth order 
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Runge-Kutta method for this system of differential equations to implement a ray tracing 
algorithm by reversing the proper time evolution. 

In order to input the metric components g a p in these equations, we have used numerical 
tables of the metric functions a, 7, p and u, generated by the code RNS with a high resolution 
mesh (1001 x 501) of the compactified quasi-isotropic radial coordinate s = f/(f + f e ), 
and the angular coordinate p = cos 8. Since the numerical scheme of the ray tracing code 
needs to evaluate the functions g a p = g a p(a(s, p),j(s, p), p(s, p),u(s, p)), defined over the 
2-dimensional grid (s,p) at any arbitrary point (f, 9), we applied bi-cubic interpolation 
between the discrete points of the numerical tables. The derivatives with respect to r and 
9 (or f and 9) in the evolution equations were replaced by finite differences, with a fine 
mesh, which we could calculate in any point from the quasi continuous interpolated function 
values. For a given set of canonical data (f, 9 = 0,pf,pe), we traced the (anti)neutrinos back 
from the rotational axis to the equatorial plane where they are supposed to be produced. 
The parameter uq was set to an arbitrary value, e.g. 0.9c, and the condition of the vanishing 
angular momentum, L = 0, gives the trajectories that intersect the rotational axis. 

With these values for the constants of motion we can determine the magnitude of p? 
from the norm p ■ p = 0, and we have a consistent set of canonical data at the axis for a 
given pq. Since Eq. f )20|) determines the relation between pg and the collision angle 9 U , the 
ray tracing code can provide the radius R at which a null particle hitting the axis at a given 
radius f and with a given collision angle left the surface of the disk in the equatorial plane. 

APPENDIX C: PHYSICAL PARAMETERS OF NEUTRON AND QUARK 
STARS 

Here we present all the data for the physical parameters of the compact stars with total mass 
M ~ 1.8M Q , and different rotational frequencies, which were obtained by the RNS code. 
In the Tables p c is the central density, M is the gravitational mass, M is the rest mass, 
R e is the circumferential radius at the equator, Q is the angular velocity, Q p is the angular 
velocity of a particle in circular orbit at the equator, T/W is the rotational-gravitational 
energy ratio, J is the angular momentum, / is the moment of inertia, $2 gives the mass 
quadrupole moment, h + is the height from the surface of the last stable co-rotating circular 
orbit in the equatorial plane, h_ is the height from surface of the last stable counter-rotating 
circular orbit in the equatorial plane, Z p is the polar redshift, Zj, is the backward equatorial 
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resdhift, Zj is the forward equatorial redshift, u c /Q is the ratio of the central value of the 
potential u to Q, r e is the coordinate equatorial radius, and r p /r e is the axes ratio (polar to 
equatorial), respectively. 

We have obtained Ri n from the formula Ri n = R e + h + , 1.6., cLS cL the sum of the stellar 
radius in the equator and the height from the surface of the last stable co-rotating circular 
orbit in the equatorial plane. Thus the inner edge of the accretion disk is located at the 
radius of the innermost stable circular orbit r ms , or at the equatorial radius R e of the stellar 
surface, if h + is zero. For the RMF stiff, STOSO, Q and CFL150 type EOS the code could 
not provide a physically reasonable values for h+, and we h ave used a truncated form of the 
analytical approximation given by ( IShibata fc Sasaki 19981 ). 

r ms « 3r fl (l- 0.54433a* - 0.22619a;; + 0.17989Q 2 - 0.23002^+ 

0.26296a*Q 2 - 0.29693a 4 + 0.44546a^Q 2 ) , (CI) 

where a* = cJ/GM 2 and Q 2 = c 4 $ 2 /G 2 M 3 are the dimensionless spin parameter and the 
dimensionless mass quadrupole moment, defined in terms of the spin angular momentum J 
and the mass quadrupole moment $ 2 . The results were checked by calculating the zeros of 
the second order derivative of the effective potential V with respect to the radial coordinate. 
For the calculation we applied a polynomial fit of the metric components in the equatorial 
plane, and computed the first and the second order derivatives of the polynomials, both 
analytically and numerically. Since the analytic forms of the derivatives of the polynomial 
fit were proved to be a good fit for the numerical derivatives too, we could use them to 
calculate V rr . The zeros of V rr were in an agreement with the results of Eq. (1C1I) within 
10%. 
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1.8043 


1.8076 


1.8022 


1.8054 


1.8071 


M [M Q ] 


2.0701 


2.0689 


2.0649 


2.0632 


2.0627 


2.0664 


2.0687 


2.0582 


2.0589 


2.0570 


-Re [km] 


11.2738 


11.2942 


11.3274 


11.3734 


11.4277 


11.5055 


11.5771 


11.7052 


11.8398 


12.0091 


J2[10 3 s _1 


] 0.5246 


1.0481 


1.5119 


2.0268 


2.5064 


3.0834 


3.5030 


4.0075 


4.5071 


4.9978 


n p [10 3 s- 


x ] 12.8662 


12.7900 


12.6968 


12.5925 


12.4873 


12.3606 


12.2508 


12.0334 


11.8445 


11.6131 


T/W[W~ 


2 ] 0.0329 


0.1317 


0.2764 


0.5016 


0.7753 


1.1888 


1.5535 


2.0921 


2.7103 


3.4355 


cJ/GM^ 


0.1109 


0.2219 


0.3208 


0.4324 


0.5385 


0.6710 


0.7707 


0.8909 


1.0198 


1.1535 


/[10 45 g cm 2 ] 1.8570 


1.8611 


1.8651 


1.8749 


1.8883 


1.9126 


1.9336 


1.9538 


1.9885 


2.0283 


4> 2 [10 42 g cm 2 ] 0.7995 


3.0606 


6.3866 


11.5847 


17.9592 


27.6936 


36.3965 


49.5047 


64.8507 


83.3478 


h+ [km] 


4.3861 


4.0805 


3.7865 


3.4791 


3.2005 


2.8939 


2.6645 


2.2884 


1.9918 


1.6754 


/i_[km] 


4.9786 


5.2653 


5.4976 


5.7775 


6.0495 


6.4120 


6.6746 


6.8938 


7.1985 


7.4826 


Zp 


0.3761 


0.3761 


0.3753 


0.3755 


0.3762 


0.3787 


0.3806 


0.3788 


0.3807 


0.3820 


z s 


0.3390 


0.3021 


0.2690 


0.2334 


0.2006 


0.1623 


0.1342 


0.0969 


0.0620 


0.0263 


Zh 


0.4137 


0.4516 


0.4847 


0.5234 


0.5607 


0.6087 


0.6446 


0.6837 


0.7287 


0.7737 


U! c /Q 


0.5987 


0.5981 


0.5965 


0.5954 


0.5946 


0.5948 


0.5947 


0.5901 


0.5888 


0.5865 


r e [km] 


8.4046 


8.4247 


8.4604 


8.5045 


8.5543 


8.6196 


8.6815 


8.8138 


8.9366 


9.0959 


r p /r e 


0.9989 


0.9955 


0.9906 


0.9830 


0.9738 


0.9600 


0.9479 


0.9300 


0.9097 


0.8860 


Table CI. The physical parameters for a neutron star model with a DH type EOS, a total mass of 1.8Mq 
from 5 X lO 2 ^ 1 to 5 X lO^" 1 . 


and angular velocities 




p c [10 15 g/cm 3 ] 


1.2640 


1.2745 


1.2710 


1.2650 


1.2500 


1.2400 


1.2200 


1.2130 






M[M Q ] 


1.7788 


1.7972 


1.7996 


1.8019 


1.8038 


1.8060 


1.8000 


1.8083 






M [M ] 


2.0746 


2.0988 


2.1006 


2.1016 


2.1013 


2.1016 


2.0895 


2.0976 






i? e [km] 


10.9336 


10.9627 


11.0009 


11.0542 


11.1299 


11.2026 


11.3140 


11.4073 






Q[10 3 s- 1 ] 


1.5271 


2.0791 


2.5175 


3.0120 


3.5665 


4.0060 


4.5288 


4.9377 






^[H^s -1 ] 


13.3079 


13.2923 


13.2144 


13.1141 


12.9808 


12.8637 


12.6656 


12.5504 






T/W[W- 2 ] 


0.2913 


0.5373 


0.7936 


1.1478 


1.6342 


2.0917 


2.7466 


3.3127 






cJ/GMl 


0.3215 


0.4459 


0.5438 


0.6564 


0.7860 


0.8928 


1.0189 


1.1312 






/[10 45 g cm 2 ] 


1.8505 


1.8847 


1.8985 


1.9152 


1.9369 


1.9586 


1.9772 


2.0134 






$ 2 [10 42 gcm 2 ] 


6.7680 


12.4757 


18.4095 


26.6997 


38.1889 


49.1323 


64.9632 


79.0479 






h+ [km] 


3.9923 


3.8447 


3.6152 


3.3529 


3.0510 


2.8164 


2.4695 


2.2982 






h- [km] 


5.7258 


6.2159 


6.4911 


6.8003 


7.1428 


7.4236 


7.6828 


8.0085 






Zp 


0.3901 


0.3980 


0.3996 


0.4014 


0.4033 


0.4052 


0.4042 


0.4086 








0.2844 


0.2532 


0.2244 


0.1918 


0.1550 


0.1258 


0.0886 


0.0623 






z b 


0.4993 


0.5495 


0.5847 


0.6253 


0.6717 


0.7102 


0.7524 


0.7944 








0.5814 


0.5876 


0.5878 


0.5878 


0.5874 


0.5870 


0.5834 


0.5848 






r e [km] 


8.0918 


8.0847 


8.1153 


8.1599 


8.2259 


8.2888 


8.4026 


8.4738 






r p /r e 


0.9907 


0.9830 


0.9750 


0.9640 


0.9490 


0.9350 


0.9150 


0.8980 





Table C2. The physical parameters for a neutron star model with APR EOS, a total mass of 1.8Mq and angular velocities 
from 1.5 X 10 3 s- 1 to 5 X H^s" 1 . 
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p c [10 15 g/cm 3 ] 


3.4000 


3.2000 


3.1000 


3.0000 


2.9000 


2.7000 


2.5000 


2.3500 


2.1800 


2.1500 


M [M ] 


1.7891 


1.7925 


1.7946 


1.7976 


1.8002 


1.8018 


1.8006 


1.8003 


1.8067 


1.8049 


M [M ] 


2.1026 


2.1073 


2.1098 


2.1130 


2.1155 


2.1159 


2.1118 


2.1083 


2.1107 


2.1077 


-R e [km] 


9.3362 


9.4389 


9.4986 


9.5670 


9.6384 


9.7808 


9.9380 


10.0843 


10.3195 


10.3498 


f7[10 3 s- 1 ] 


0.5525 


1.0536 


1.5232 


2.0716 


2.5101 


3.0357 


3.5133 


4.0414 


4.9989 


5.0398 


£} p [10 3 s _1 ] 


17.0064 


16.6910 


16.4982 


16.2884 


16.0871 


15.7139 


15.3184 


14.9718 


14.4799 


14.4116 


T/W[10-' 2 } 


0.0202 


0.0764 


0.1633 


0.3095 


0.4659 


0.7178 


1.0182 


1.4153 


2.3168 


2.3808 


cJ/GM^ 


0.0873 


0.1702 


0.2490 


0.3435 


0.4221 


0.5237 


0.6215 


0.7317 


0.9427 


0.9536 


/[10 45 g cm 2 ] 


1.3887 


1.4195 


1.4368 


1.4572 


1.4781 


1.5160 


1.5547 


1.5912 


1.6573 


1.6629 


$ 2 [10 42 gcm 2 ] 


0.3292 


1.1571 


2.4929 


4.7963 


7.3192 


11.6367 


17.0855 


24.4570 


41.7102 


43.1197 


h+ [km] 


6.2663 


5.9751 


5.7298 


5.4476 


5.2036 


4.8317 


4.4388 


4.0492 


3.4490 


3.3853 


/i_[km] 


6.7367 


6.8894 


7.0655 


7.2848 


7.4555 


7.6162 


7.7352 


7.9136 


8.3595 


8.3522 




0.5173 


0.5094 


0.5060 


0.5032 


0.5002 


0.4916 


0.4815 


0.4748 


0.4718 


0.4695 


z s 


0.4779 


0.4348 


0.3986 


0.3578 


0.3245 


0.2801 


0.2379 


0.1954 


0.1257 


0.1208 


Zh 


0.5571 


0.5859 


0.6173 


0.6559 


0.6866 


0.7186 


0.7456 


0.7812 


0.8597 


0.8604 


U! c /Q 


0.7428 


0.7328 


0.7276 


0.7226 


0.7172 


0.7050 


0.6912 


0.6805 


0.6694 


0.6666 


r e [km] 


6.4255 


6.5254 


6.5820 


6.6451 


6.7119 


6.8525 


7.0128 


7.1580 


7.3745 


7.4081 


r p /r e 


0.9993 


0.9975 


0.9947 


0.9900 


0.9850 


0.9770 


0.9675 


0.9550 


0.9270 


0.9250 



Table C3. The physical parameters for a neutron star model with BBBAV14 EOS, a total mass of 1.8Mq and angular velocities 
from 5 X lO 2 ^ 1 to 5 X lO^" 1 . 



p c [10 15 g/cm 3 ] 


1.9000 


1.8900 


1.8800 


1.8700 


1.8500 


1.8300 


1.8000 


1.7500 


1.7000 


M [M ] 


1.8017 


1.8016 


1.8033 


1.8049 


1.8055 


1.8072 


1.8085 


1.8014 


1.7989 


M [M ] 


2.1170 


2.1161 


2.1171 


2.1180 


2.1169 


2.1170 


2.1157 


2.1027 


2.0949 


Re [km] 


10.1829 


10.2058 


10.2378 


10.2697 


10.3218 


10.3809 


10.4664 


10.5682 


10.7068 


n[10 3 s- 1 ] 


1.0141 


1.5161 


2.0785 


2.5091 


3.0185 


3.5098 


4.0696 


4.4813 


5.0407 


n p [10 3 s- 1 ] 


14.9378 


14.8458 


14.7434 


14.6556 


14.5239 


14.3900 


14.2075 


13.9766 


13.7052 


T/W[W- 2 ] 


0.0957 


0.2150 


0.4071 


0.5977 


0.8768 


1.2028 


1.6522 


2.0679 


2.7144 


cJ/GMl 


0.1897 


0.2842 


0.3918 


0.4757 


0.5768 


0.6771 


0.7952 


0.8832 


1.0106 


/[10 45 g cm 2 ] 


1.6439 


1.6473 


1.6567 


1.6662 


1.6792 


1.6954 


1.7173 


1.7322 


1.7620 


<J> 2 [10 42 gcm 2 ] 


1.8153 


4.0693 


7.6582 


11.2359 


16.5747 


22.8810 


31.7467 


40.2724 


53.8111 


/)+[km] 


5.2697 


5.0071 


4.7255 


4.5082 


4.2300 


3.9658 


3.6472 


3.3156 


2.9369 


h- [km] 


6.2820 


6.5233 


6.8101 


7.0316 


7.2789 


7.5282 


7.8059 


7.9256 


8.1701 


Zp 


0.4478 


0.4477 


0.4486 


0.4495 


0.4497 


0.4505 


0.4509 


0.4461 


0.4441 


Zf 


0.3764 


0.3415 


0.3036 


0.2747 


0.2399 


0.2068 


0.1686 


0.1363 


0.0956 


z b 


0.5209 


0.5576 


0.6007 


0.6346 


0.6744 


0.7144 


0.7605 


0.7885 


0.8337 




0.6546 


0.6536 


0.6531 


0.6526 


0.6509 


0.6494 


0.6469 


0.6400 


0.6343 


r e [km] 


7.2806 


7.3023 


7.3290 


7.3559 


7.4037 


7.4561 


7.5343 


7.6451 


7.7814 


r p /r e 


0.9969 


0.9931 


0.9870 


0.9810 


0.9722 


0.9620 


0.9480 


0.9350 


0.9150 



Table C4. The physical parameters for a neutron star model with BBBParis EOS with a total mass of 1.8M and angular 
velocities from 10 3 s _1 to 5 X 10 3 s _1 . 
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Pc [10 15 g/cm 3 ] 


0.7290 


0.7250 


0.7200 


0.7100 


0.7000 


0.6850 


0.6670 


0.6400 


0.6100 


0.5700 




M[M Q ] 


1.8065 


1.8053 


1.8069 


1.8062 


1.8078 


1.8054 


1.8062 


1.8000 


1.7993 


1.8153 




M [M ] 


2.0443 


2.0421 


2.0429 


2.0399 


2.0395 


2.0336 


2.0304 


2.0172 


2.0092 


2.0173 




Re [km] 


13.5283 


13.5644 


13.6313 


13.7402 


13.8690 


14.0370 


14.2746 


14.6091 


15.0696 


15.9090 




^[io^- 1 ] 


0.5335 


0.9922 


1.5077 


2.0600 


2.5368 


3.0031 


3.5078 


4.0159 


4.5104 


5.0740 




n P [io 3 s- 1 ] 


9.7975 


9.7329 


9.6510 


9.5304 


9.4096 


9.2539 


9.0579 


8.7849 


8.4531 


7.9384 




T/W[W- 2 ] 


0.0679 


0.2368 


0.5533 


1.0555 


1.6387 


2.3739 


3.3823 


4.7272 


6.4613 


9.2817 




cJ/GMl 


0.1641 


0.3065 


0.4702 


0.6512 


0.8160 


0.9849 


1.1851 


1.4065 


1.6653 


2.0742 




/[10 45 g cm 2 ] 


2.7037 


2.7147 


2.7411 


2.7782 


2.8268 


2.8824 


2.9691 


3.0780 


3.2448 


3.5927 




<J> 2 [10 42 gcm 2 ] 


2.8764 


10.1367 


23.4709 


44.8791 


70.2820 


102.8550 


148.8100 


212.4910 


299.7310 


458.3680 




h+ [km] 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 




h_ [km] 


0.0000 


0.0000 


0.0000 


— oo 


— oo 


— oo 


— oo 


— oo 


— oo 


— oo 




Zp 


0.2852 


0.2852 


0.2861 


0.2867 


0.2880 


0.2884 


0.2898 


0.2896 


0.2905 


0.2945 






0.2456 


0.2118 


0.1746 


0.1342 


0.0994 


0.0641 


0.0252 


-0.0170 


-0.0609 


-0.1178 




z b 


0.3251 


0.3600 


0.4008 


0.4453 


0.4859 


0.5260 


0.5726 


0.6205 


0.6734 


0.7492 




ui c /Q 


0.4748 


0.4741 


0.4740 


0.4727 


0.4720 


0.4697 


0.4676 


0.4628 


0.4584 


0.4549 




r e [km] 


10.6958 


10.7316 


10.7913 


10.8947 


11.0128 


11.1748 


11.3977 


11.7252 


12.1652 


12.9435 




r p /r e 


0.9977 


0.9920 


0.9815 


0.9650 


0.9460 


0.9225 


0.8910 


0.8500 


0.7990 


0.7200 




Table C5. The physical parameters for a neutron star model with RMF stiff EOS, a total mass of 1.8Mq and 
from 5 X 10 2 s- 1 to 5 X lO^" 1 . 


angular velocities 




p c [10 15 g/cm 3 ] 


0.4800 


0.4770 


3.2000 


0.4750 


0.4700 


0.4620 


0.4540 


0.4400 


0.4230 


0.3980 


0.3690 


M[M Q ] 


1.8061 


1.8031 


1.7925 


1.8076 


1.8086 


1.8017 


1.8088 


1.8002 


1.8039 


1.8033 


1.8541 


Mo [M ] 


2.0118 


2.0070 


2.1073 


2.0113 


2.0103 


1.9990 


2.0034 


1.9876 


1.9842 


1.9725 


2.0155 


Re [km] 


15.2003 


15.2651 


9.4389 


15.3585 


15.5171 


15.7140 


16.0281 


16.4469 


17.1311 


18.2747 


21.0394 


n[10 3 s- 1 ] 


0.5081 


1.0560 


1.0536 


1.5041 


2.0308 


2.5000 


3.0348 


3.5301 


4.0582 


4.5293 


4.9020 


n p [10 3 s- 1 ] 


8.2289 


8.1525 


16.6910 


8.0845 


7.9708 


7.8258 


7.6457 


7.3877 


7.0275 


6.4676 


5.3725 


T/W[W- 2 ] 


0.0865 


0.3780 


0.0764 


0.7757 


1.4473 


2.2664 


3.4814 


5.0168 


7.2439 


10.3451 


14.9104 


cJ/GMl 


0.1891 


0.3949 


0.1702 


0.5701 


0.7835 


0.9799 


1.2345 


1.4877 


1.8274 


2.2429 


2.9524 


/[10 45 g cm 2 ] 


3.2706 


3.2868 


1.4195 


3.3309 


3.3905 


3.4448 


3.5750 


3.7037 


3.9573 


4.3520 


5.2932 


$ 2 [10 42 gcm 2 ] 


4.8687 


21.1326 


1.1571 


43.4479 


81.4042 


127.9780 


199.1870 


290.8950 


432.8830 


650.3940 


1061.9200 


ft 4- [km] 


— oo 


0.0000 


5.9751 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


h- [km] 


1.3339 


1.9240 


6.8894 


2.4997 


3.1777 


3.7552 


4.5759 


5.2774 


6.2012 


7.1106 


7.8501 


Zp 


0.2415 


0.2416 


0.5094 


0.2432 


0.2446 


0.2448 


0.2480 


0.2486 


0.2515 


0.2529 


0.2601 




0.2020 


0.1596 


0.4348 


0.1260 


0.0856 


0.0481 


0.0056 


-0.0380 


-0.0885 


-0.1461 


-0.2296 


z b 


0.2814 


0.3251 


0.5859 


0.3634 


0.4094 


0.4504 


0.5038 


0.5534 


0.6165 


0.6842 


0.7891 


u) c /Q 


0.4201 


0.4194 


0.7328 


0.4203 


0.4202 


0.4183 


0.4191 


0.4160 


0.4145 


0.4099 


0.4102 


r e [km] 


12.3910 


12.4567 


6.5254 


12.5372 


12.6848 


12.8822 


13.1688 


13.5832 


14.2357 


15.3506 


18.0144 


r-p/r-e 


0.9970 


0.9870 


0.9975 


0.9735 


0.9510 


0.9240 


0.8850 


0.8370 


0.7700 


0.6800 


0.5400 



Table C6. The physical parameters for a neutron star model with STOS0 EOS, a total mass of 1.8Mq and angular velocities 
from 5 x 10 2 S - 1 to 5 x 10 3 s- 1 . 
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p c [10 15 g/cm 3 ] 


1.0800 


1.0700 


1.0500 


1.0300 


1.0000 


0.9800 


0.9490 


0.9310 


M [M ] 


1.8039 


1.8065 


1.8032 


1.8029 


1.8000 


1.8014 


1.7999 


1.8071 


M [M Q ] 


2.1221 


2.1241 


2.1176 


2.1146 


2.1071 


2.1057 


2.0989 


2.1045 


Re [km] 


11.2268 


11.2717 


11.3243 


11.3914 


11.4900 


11.5804 


11.7097 


11.8329 


^[10 3 s- 1 ] 


1.4909 


2.0426 


2.5022 


2.9918 


3.5684 


3.9966 


4.5072 


4.9097 


n P [io 3 s- 1 ] 


12.8593 


12.7703 


12.6630 


12.5507 


12.3985 


12.2828 


12.1180 


11.9932 


T/W[W- 2 ] 


0.3678 


0.6975 


1.0654 


1.5517 


2.2722 


2.9149 


3.8377 


4.6703 


cJ/GMl 


0.3650 


0.5047 


0.6224 


0.7523 


0.9103 


1.0352 


1.1904 


1.3280 


/[10 45 g cm 2 ] 


2.1519 


2.1717 


2.1862 


2.2100 


2.2420 


2.2765 


2.3212 


2.3772 


$ 2 [10 42 gcm 2 ] 


10.3179 


19.1974 


30.1157 


43.9566 


64.5770 


83.4277 


110.8010 


136.1140 


h+ [km] 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


h_ [km] 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


Zp 


0.3822 


0.3837 


0.3832 


0.3839 


0.3843 


0.3859 


0.3870 


0.3906 




0.2775 


0.2405 


0.2084 


0.1752 


0.1355 


0.1066 


0.0712 


0.0443 


z b 


0.4905 


0.5339 


0.5685 


0.6080 


0.6550 


0.6930 


0.7386 


0.7805 


U> c /Q 


0.5400 


0.5399 


0.5374 


0.5359 


0.5329 


0.5317 


0.5289 


0.5292 


r e [km] 


8.3558 


8.3912 


8.4438 


8.5041 


8.5967 


8.6746 


8.7925 


8.8900 


r p /r e 


0.9893 


0.9800 


0.9690 


0.9549 


0.9350 


0.9170 


0.8920 


0.8700 



Table C7. The physical parameters for a quark star model with Q EOS, a total mass of 1.8Mq and angular velocities from 
1.5 X 10 3 s- 1 to 5 X 10 3 s- 1 . 



p c [10 15 g/cm 3 ] 


0.7830 


0.7790 


0.7750 


0.7670 


0.7580 


0.7470 


0.7350 


0.7200 


0.7050 


M [M ] 


1.8003 


1.8000 


1.8036 


1.8015 


1.8012 


1.7994 


1.8022 


1.7997 


1.8069 


M [M e ] 


2.1162 


2.1149 


2.1179 


2.1129 


2.1098 


2.1042 


2.1035 


2.0954 


2.0986 


fi e [km] 


11.7280 


11.7561 


11.8075 


11.8591 


11.9321 


12.0175 


12.1406 


12.2707 


12.4636 


^[io^- 1 ] 


1.0163 


1.4994 


2.0421 


2.5147 


3.0197 


3.5035 


4.0415 


4.5235 


5.0640 


Q p [10 3 s- 1 ] 


12.0724 


12.0028 


11.9281 


11.8403 


11.7438 


11.6396 


11.5143 


11.3772 


11.2109 


T/W[W- 2 ] 


0.1987 


0.4351 


0.8147 


1.2517 


1.8350 


2.5197 


3.4382 


4.4366 


5.7661 


cJ/GMl 


0.2749 


0.4071 


0.5603 


0.6944 


0.8430 


0.9894 


1.1645 


1.3263 


1.5338 


/[10 45 g cm 2 ] 


2.3770 


2.3865 


2.4112 


2.4269 


2.4533 


2.4820 


2.5323 


2.5767 


2.6619 


* 2 [10 42 gcm 2 ] 


7.1160 


14.9170 


27.8321 


42.6462 


62.2833 


85.5102 


117.0170 


151.5570 


198.8980 


/i+ [km] 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


h- [km] 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


Zp 


0.3539 


0.3545 


0.3566 


0.3573 


0.3589 


0.3604 


0.3635 


0.3654 


0.3703 




0.2819 


0.2487 


0.2125 


0.1803 


0.1463 


0.1135 


0.0775 


0.0441 


0.0070 


z b 


0.4276 


0.4641 


0.5077 


0.5451 


0.5872 


0.6285 


0.6784 


0.7234 


0.7813 




0.5010 


0.5008 


0.5016 


0.5008 


0.5005 


0.4997 


0.4999 


0.4988 


0.4997 


r e [km] 


8.8688 


8.8932 


8.9318 


8.9796 


9.0431 


9.1197 


9.2223 


9.3399 


9.4977 


r p /r e 


0.9935 


0.9864 


0.9750 


0.9620 


0.9450 


0.9250 


0.8990 


0.8710 


0.8350 



Table C8. The physical parameters for a neutron star model with CFL150 EOS, a total mass of 1.8Mq and angular velocities 
from 10 3 s" 1 to 5 X 10 3 s- 1 . 
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